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1. Introduction 

Lattice QCD opened the possibility of providing a model independent ab initio calculation 
of the QCD mass spectrum. Extracting signals from excited states has however proven to be a 
formidable task. A previous study of the BGR-Collaboration [|I|] employed the variational tech- 
nique [§,|3|] to study excited mesons on the lattice using standard meson interpolators with Gaussian- 
smeared sources and sinks. 

For an accurate representation of excited states in the variational approach it is crucial to use a basis 
with good overlap with the physical states. We therefore construct additional meson-interpolators 
with derivative sources obtained by a covariant derivative acting on the Jacobi smeared sources and 
explore their effect in the variational approach. 

2. Details of the calculation 

2.1 The variational method 

We use a basis of interpolators Oi,i = l,...,N with the quantum numbers of the desired states 
(projected to zero momentum) and compute the matrix of cross correlations: 

C(t)ij = (Oi(t)O](0)) =^(0\Oi\n)(n\o)\0)e- tM " . (2.1) 

n 

It can be shown [|3j] that the solutions to the generalized eigenvalue problem 

C{t)vi = Xi(t)C(t )vi , (2.2) 

behave as 

Xi(t) oc e - ,M ' (1 + (e- tAM >)) , (2.3) 

where AM, is the mass difference between the state i and the closest lying state. 
Therefore, up to the given order, the largest eigenvalue decays with the mass of the ground state, 
the second largest eigenvalue with the mass of the first excited state and so on. For an accurate 
description, the interpolators should be linearly independent, as orthogonal as possible and possess 
a strong overlap with physical states. The task will be to construct new interpolators which are 
independent of the standard interpolators and have a good overlap with excited states. 

2.2 Effective masses and eigenvectors 

To display the results, one commonly plots the effective masses determined from ratios of 
eigenvalues 

"«H) =h (>li))' <24) 

A fit to the exponential decay of the eigenvalues is then performed in the time interval where the 
effective mass plot shows a plateau. As a consistency check, the entries of the corresponding eigen- 
vectors should also show plateaus and can serve as a fingerprint for the "wave functions". To deal 
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with eigenvectors that are orthonormal we calculate the eigenvectors of the modified eigenvector 
problem 

C(/ )-*C(0C(U))-^ = WW ■ (2.5) 
(This assumes C(to) to be positive definite.) 

2.3 Smeared sources and sinks 

The first step for all our sources is Jacobi smearing [Q, |5|] of point sources so located at timeslice 
t = 0: 

s { o a ' a \y,0) P ,c = S(y,0)8 pa 8 ca , (2.6) 

N 

s (a,a) = £ ^ H n ? (2 7) 

n=0 
3 

H{x,y) = £ {Ui{x,0)d{x + ly) + Ui{x-i^8{x-ly)) . (2.8) 



i=i 



In two different parameter values for k and N, giving rise to wide and narrow sources (S w ,S n ), 
have been used. In addition, we now construct covariant derivatives which act upon a wide smeared 
source to form our derivative quark sources W c /, : 

Pi(x,y) = Ui(x,0)8{x + i,y) - Ui(x-t,0) f 8(x-i,y) , (2.9) 
W d .=PiS w . (2.10) 

With these sources, meson interpolators of definite quantum numbers are constructed. 

2.4 Interpolators used 

Table 1 shows the interpolators used for different meson channels. The ones labeled "old 
interpolators" are those of [|I|], while the "new interpolators" contain at least one derivative. The 
total number of interpolators as given in the last column can be obtained by combining all non- 
degenerate combinations of the different smearings with all possible Dirac structures. 

In some cases, an (anti-) symmetrization of the interpolators is necessary to obtain the cor- 
rect behavior under charge conjugation. Therefore, interpolators denoted as u c iTd n / w in the table 
should be read as U c iTd n / w — u n / w Tdd r We restrict ourselves to light, isovector (1=1) mesons with 
degenerate quark masses m u = m^. 

2.5 Technicalities 

For our analysis we used 99 uncorrected quenched gauge configurations with Chirally Im- 
proved (CI) H ^] quarks and the Liischer-Weisz gauge action [Q]. We work on a 16 3 x 32 lattice 



with a = 0.148 fm determined from [1C] the Sommer parameter. We study several quark mass pa- 
rameters in the range am q = 0.02 . . . 0.2. All errors we quote are statistical errors determined with 
the jackknife method. 
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Table 1: Meson interpolators; n/w denote narrow and wide Jacobi smearing and di stands for a derivative 
source in the /-direction. The last column shows the total number of different interpolators. 

3. Results 

To get an idea about the properties of the correlators, we first plot the normalized diagonal 
entries of the correlation matrices for the pseudoscalar and vector mesons. In figure 1 the simple 
interpolators have been colored blue, while the new correlators from derivative sources appear in 
red. As can be seen from the behavior at small euclidean times t, the new correlators have stronger 
contributions from excited states, i.e., they start out with a steeper slope. However, all correlators 
are dominated by the ground state at large time separations. 

Figure 2 shows a comparison of selected effective mass plateaus from the optimal set of old 
interpolators compared to the new results (obtained from the optimal combination of both old and 
new interpolators) at different values of the quark mass. Plateaus are also found in the components 
of the corresponding eigenvectors. 

For the pseudoscalar mesons (graphs on the upper l.h.s.) we obtain improved results for the 
excited states at all quark masses. Both the old and the new sets are obtained with a basis of four 
interpolators. This shows that the observed effect does not just stem from simply enlarging the 
basis. Some interpolators mainly contribute noise without significantly enhancing the overlap to 
the physical states. Looking at the modified eigenvectors of equation (2.5) as depicted in figure 3, 
clear plateaus are observed and the new interpolators (second and third entries) mainly contribute 
to the second excited state which could not be observed using the old set of interpolators alone. 

Similarily for the vector mesons (upper r.h.s. of figure 2), improved results are obtained for 
the excited states at all quark masses. Here, however, the new basis consists of six interpolators 
while the old basis only contains four. Both the ground state and newly observed third excited 
state are dominated by the old interpolators while the first and second excited states are dominated 
by the new interpolators.We also would like to note that the finer (a 0.12fm) quenched lattice 
considered in |J|] allowed for a cleaner separation of the first and second excited states. 

For the scalar (lower l.h.s. of figure 2) and pseudovector (lower r.h.s. of figure 2) mesons, there 
is a noticeable improvement of the ground state plateaus. The old graph for the pseudovectors 
stems from a basis of three interpolators while the new graph is obtained from just one correlator 
alone. The location of the plateau agrees within error bounds but the errors obtained from the new 
interpolator are considerably smaller. For the scalars, the plots have been obtained from a basis 
of three interpolators. Using derivative sources, plateaus can be extracted at small values of the 
quark mass where no plateau is observed in the old interpolators. For the quenched data, scalar 
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Figure 1: Diagonal correlators for the pseudoscalar (top) and vector (bottom) mesons at am q = 0.04. n/w 
denote narrow and wide Jacobi smearing and di stands for a derivative source in the /-direction. 

correlators contain contaminations from ghost states and it has been shown, that the variational 
method can isolate such unphysical contributions [Q]. Overall, ghost contributions to derivative 
interpolators seem to be smaller and some interpolators show almost no sign of such contributions. 
Where present, the variational method separates the ghosts from the physical signal and the ground 
state plateau can be clearly extracted (shown in figure 2). Concerning a full dynamical simulation, 
data from only a few dynamical configurations suggests that interpolators of the type Ud.yid n / W , 
which are important in the quenched case, seem to have very limited overlap in the dynamical case. 

4. Conclusions and upcoming investigations 

In this article we have presented first results for pseudoscalar, scalar, vector and pseudovector 
mesons with derivative quark sources and CI fermions on a quenched lattice. We employed the 
variational technique with both Jacobi-smeared sources of a Gaussian type as well as derivative- 
type sources. 

It has been demonstrated that interpolators constructed with derivative quark sources lead to an 
enhanced signal for a variety of different meson channels. For pseudoscalar and vector mesons, the 
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Figure 2: Effective mass plots for the pseudoscalar (upper left), vector (upper right), scalar (lower left) and 
pseudovector (lower right) mesons. We compare the results from the best selection of simple interpolators 
("old") to the case where the new derivative interpolators are added ("new"). The horizontal lines indicate 
the fit ranges and their error bars indicate the statistical error of the fit. 
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Figure 3: Eigenvector components of the modified eigenvector problem for the pion ground state (l.h.s. 
plot), first and second excited states (center and r.h.s. plot). The plots correspond to the effective mass 
plateaus shown in figure 2. 

correlators constructed from such interpolators display a significantly better overlap with excited 
states. In both cases excited states could be observed that were not found with the smaller basis. 
For the scalar and pseudovector mesons, overlap with the ground states can be improved which 
leads to ground state masses with significantly smaller statistical errors. 

We have clearly demonstrated that derivative quark sources provide an essential enlargement of 
the basis of interpolators in excited meson spectroscopy. Although we are convinced that derivative 
interpolators will be an important tool also in the dynamical case, it is probably too early to claim 
a one to one correspondence of the interpolator content. 

We are currently extending our calculations to dynamical configurations obtained in the BGR 
dynamical CI project. Furthermore, we investigate the construction of baryon interpolators from 
derivative quark sources and the effects of link smearing on our interpolators. 
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